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Abstract 

We obtain the exact solution for the Burgers equation with a time depen- 
dent forcing, which depends linearly on the spatial coordinate. For the case 
of a stochastic time dependence an exact expression for the joint probability 
distribution for the velocity fields at multiple spatial points is obtained. A 
connection with stretched vortices in hydrodynamic flows is discussed. 

PACS number(s): 02.50.-r, 05.40.-a, 05.90.-a, 47.32. Cc 

Burgers equation 0, [2], [H] has extensively been investigated as a model 
equation for the study of the problem of turbulence. The homogeneous equa- 
tion can be exactly solved by a Hopf-Cole transformation [3], [3] such that 
many properties of the spatio-temporal behaviour can be obtained by an- 
alytical means. In recent years, the stochastically forced Burgers equation 
has attracted much interest due to the work of Polyakov [Oj. He considered 
a white noise forcing f(x,t) with spatial correlations < f(x,t)f(x',t) >= 
)8(t — t'). This equation has been found to exhibit scaling behaviour. 
Polyakov suggested to use an operator product expansion in order to assess 
the scaling behaviour of the velocity increments. 

The purpose of the present note is to present the general solution of the 
forced Burgers equation 



d d d 2 

L - + ,M)-Mr,t) = v - 

where G(t) is an arbitrary function of time t. As we shall discuss below, 



l^ + v(r,t)-]v(r,t) = u—v(r,t) + G(t)r , (1) 



Vahhart@uni-muenster.de, Tel.:+49-251-8334924, Fax: +49-251-8336328 



1 



this problem is related to the Polyakov problem. The investigation of the 
special type of forcing, f(r,t) = rG(t), has the advantage that one can find 
an analytic solution for an arbitrary function G(t). 

This solution can be found by direct calculation using the ansatz 

v(r,t) = a{t)r + A(t)w(A(t)r,r(t)) . (2) 

Thereby, the functions a(t), A(t), r(t) are determined by the following solv- 
ability conditions: 

a(t) + a(t) 2 = G{t) 

A = - a (t)A , 

i = A 2 . (3) 
The quantities A(t) and r(t) can be determined as a function of a(t), 

A(t) = e-fo , r = f'dt'e- 2 ^^") (4) 

Jo 

In the following we shall see that the quantity 

W - ^ = (5) 

will play a major role. 

In order that the ansatz eq. (j2J) obeys the forced Burgers equation the 
function u>(£, r) has to obey the homogeneous equation. For instance, r) 
could be the single shock solution 

w(£,t) = -A tank ^ (6) 

Of interest are also scaling solutions of the form jH] 



which e.g. evolve from the initial condition u>(£, 0) = 

Since the solution of the homogeneous Burgers equation is known in terms 
of the initial conditions the function A(t)w(£, r) can be represented according 
to 

A(t)MAr,r) = -2vA(t)^-ln f°° dr ' e Hr-r')V{^)) e - r!^m^,o)/2 V 

or J-oo 
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We mention that the structure of the solution, eq. (J2J), is similar to the 
solution describing stretched vortices in the Navier-Stokes equation (see e.g. 
0). The analogy becomes even more clear since, as we shall briefly mention 
below, the treatment can be extended to the case of the multidimensional 
Burgers equation. Stretched vortex solutions of the Navier-Stokes equation 
are supposed to play a fundamental role in the organization of turbulence. 
Stretching has the effect to render a decaying vortex like the Lamb-Oseen 
vortex into a stationary vortex (e.g. the Burgers vortex). Similar phenomena 
are expected in the present case of the forced Burgers equation. 

As can be seen from the solvability conditions, the crucial point is to 
obtain a solution of the nonlinear differential equation for a(t). The behaviour 
of a(t) can be visualized by introducing the potential V(a) = |a 3 : The 
variable a(t) behaves as an overdamped particle moving in the potential 
V(a) under the influence of the time dependent force G(t). By introducing 
the variable b(t) according to 

°« = W) (9) 

the equation of motion for a(t) can be transformed into a linear equation of 
second order (Hill's equation) 

bit) = G{t)b(t) . (10) 

In terms of b(t), the variables A(t) and r(t) are determined by 

6(0) 



A(f) 



b(t) 



1 „&(0), 2 



T ® = Jo dt 'W (U) 

In the following, we shall briefly consider some cases of functions G(t). 

The first case is a stationary forcing, G(t) = g. If we take g to be a 
positive number, one obtains two solutions for a 

_ i 

a = ±Jg~ A = e T v^ T = . (12) 

Let us first consider the case of a positive value of a = ^fg. This solution is 
linearly stable. The contribution of a shock 

- AA{t)tanh^^-r (13) 
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decays to zero. A scale invariant solution also decays, since 



and 



A(( > W(^Lr) (14) 



Aft) e -2y/m 
lim -±J= = 2y/g- = (15) 



The case of a negative value of a = —^[g is linearly unstable. Further- 
more, shocklike contributions increase in time connected with a steepening 
of the shock. Scale invariant contributions in the long time limit become 
stationary since 

A(t) e 2 ^ 



lim -^L = IJg"^^ = JlJg (16) 

such that the solution takes the form 

v(r,t) = -^r + ^2~^W(^2~^r) . (17) 

If we take g = —uj 2 to be negative, we obtain an oscillating solution for 
b(t) and a finite time singularity for a(t): 

b(t) = A cos ujt a(t) = —uj tan u>t (18) 

Furthermore, the quantity a(t) reads: 



AW / 2, 



a (t) \l T {t) ' sinlujt 



We shall now consider the case where G(t) is a stochastic white noise 
force. The Fokker-Planck equation for the probability distribution f(a, t) =< 
5 (a — a(t)) > reads 

A stationary solution to this equation can be obtained: 

fo(a) = Ne~ v ^' Q (21) 

However, since the potential V(a) is not bounded from below, this stationary 
solution does not define a probability distribution, i.e. a stationary state does 
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not exist. Due to the quadratic nonlinearity in the Langevin equation a kind 
of stochastic finite time singularity exists, where a(t) tends to a(t) — > — oo. 
This can be seen from the fact that the variable a moves like an overdamped 
particle in the unbounded potential under the influence of a stochastic force. 
As a consequence the variable A(t) increases, too. 

Although there is a finite time singularity it is quite interesting to in- 
vestigate the structure of the probability distribution of the velocity field 
v(r,t). To this end we introduce the distribution P(a, A, r;t), which obeys 
the Fokker-Planck equation 
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d „ , d <9 , 9 d n d' - , 

oi p = < m aA ~ a7 A + Ta a + w )p (22) 

Using P(a,A,T',t) the probability distribution of the velocity field v(r,t), 
f(v, r, t) =< 5(v — v(r, t) > takes the following form 

f(v,r,t) = JdAdrdaP(a,A,T,t)S(v-ar-Aw(Ar,T)) 

= l -JdrdAP C~ kW ^ r \ Kr,t) (23) 

The joint pdf for the velocity field at multiple scales r\ takes the form 

f({vi, n}, t) = f dAdrdaP(a, A, r, t)U^ =1 5( Vi - an - Aw(Ar i} r)) . (24) 

This formula demonstrates how the spatial coherence of the stochastic force 
governs the statistics of the multiple scales pdf of the velocity field. 

Using standard methods one may easily derive an evolution equation for 
the pdf f(v, r, t): 

J^(^) + v^-f(v : r,t) + ^-J^dv'f(v',r : t) 

= -v J dv' j dr'S(r - r')v'-^f(v',r';v,r;t) 

' dv 2 

This equation is the first equation in a whole chain of equations for the joint 
pdf of the velocity field at multiple scales, 

^f(vi,r h t) + J2v i -^-f(v i ,r h t) + J2-^- J^dv'J(v' i: r h t) 

= -J> J dv' J dr'5(r - r ')v'—f{v', r'; v i: n; t) 

i 

+ J2Q(^r 3 )^f(v t ,r t ,t) (26) 



+ Q(r,r)—f(v,r,t) . (25) 
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Here, the matrix Q{r^rj) is given by 

Q(n,rj) = QriTj (27) 

We present this hierarchy, since Polyakov [S| considered the same hierarchy 
for the case of a translational invariant white noise forcing. Here, we could 
obtain a solution to this chain of equations by a direct calculation in terms 
of equation (|2~^) . 

Let us now discuss the relationship of the above problem of forced Burgers 
equation to the problem considered by Polyakov jE] of translational invariant 
white noise forcing: 



d d d 2 

{— + u(x,t)—]u(x,t) = u—u(x,t) + F(x,t) 

< F(x,t)F(x',t) > = q(x - x')6(t - 1') (28) 



Of interest are the statistical properties of the velocity increments 

v x (r, t) — u(x + r, t) — u(x, t) . (29) 
The increment v x (r,t) obeys the evolution equation 

d d d 

{— + u(x,t)—]v x (r,t)+v x (r,t)—v x (r,t) 

d 2 d 2 
= u \-g-2 v *( r > t ) ~ ^3 U *( r >*)lr=o] + fx(r,t) 

f x (r,t) = F(x + r,t)-F(x,t) 

< U(r,t)f x (r',t)> = q(r-r')-q(r)-q(r') + q(0) (30) 

The additional dissipative term arises since one has the condition 

v x {r = 0,t) = . (31) 

In a comoving coordinate system, i.e. considering the increments v x M(r,t), 
where the path x(t) is defined by 

x(t)=u(x(t),t) (32) 

we obtain the forced Burgers equation for the increment 

r d d 

i-Ql + v{r,t)—}v{x,t) 

= V^ 1 "^' 11 ^^^ 5 ' (33) 
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Here, G(r,t) is given by 



G(r, t) = F{x(t) +r,t)- F{x(t),t) (34) 

Polyakov investigated the statistics of increments in the limit of small 
values of r and approximated the function g(£) by 

9(O=<?(0)-f£ 2 + - (35) 

This relationship can only hold for |£| << yJq(0)/Q, since for large values of 
£ q{Q ~~ * 0- This approximation, however, is equivalent to retain only the 
term linear in r in the expression of the fluctiuating force: 

G(r,t) = rG{t) + h.o. (36) 

The function G(t) is related to the stochastic force F(x,t) via the rela- 
tionship 

G(t) = ^F(x,t)\ x=x(t) (37) 

In this approximation we recover our starting equation, eq. (J33|) . but now 
formulated for increments. That means, that the solutions are given by a 
similar expression compared to eq. (0) 

v(r, t) = a(t)r + A{t)w{A{t)r, r(t)) (38) 

provided we take w(£, r) as increment solutions, i.e. w(£, r) = w(x + r) — 
w(£, r), where w(£, r) is a solution of the homogeneous Burgers equation. 

Let us now briefly consider the treatment of the velocity increment statis- 
tics in the translationally invariant white noise forcing case, which has been 
performed by Polyakov on the basis of a hierarchy similar to eq. ([26)1 . The 
only difference in this hierarchy is the form of the matrix Q(r», rj). It is 
defined according to 

Q(r h rj) = q(ri) + qfa) - q(0) - q{n - rj) (39) 

Taking the approximation q(r) = q(0) — Q/2r 2 the same hierarchy arises 
since in this approximation Q(ri,rj) = QriVj. Under this approximation 
the hierarchy of Polyakov is solved by the probability distributions given in 
eq. (|24|). However, as we have indicated above, the limit t — > oo of the 
solution does not exist. Therefore, the evaluation of the stationary pdf's for 
the Polyakov problem has to go beyond the approximation of the force eq. 
(|nHJ). This indicates that all calculations which have been performed so far 
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for the translationally invariant white noise forced Burgers equation exclu- 
sively based on the approximation should lead to the above probability 
distributions, which, as we have indicated above, have no stationary time 
limit. 

Let us now finally consider the multidimensional forced Burgers equation: 
|U(r, t) + v(r, t) ■ V r v(r, *) = z/Av(r, *) + G(t)r (40) 

Here, G(t) denotes a time dependent matrix. We look for solutions of the 
form 

v = A(t)r + U(r, t) (41) 

The matrix A(t) is determined in such a way as to cancel the force field. This 
leads to the matrix equation 

A + AA = G(t) , (42) 

whereas the additional field U(r, t) has to fulfill the equation 

-^U(r, t) + A ■ U(r, t) + (A ■ r) • VU(r, t) + U(r, t) ■ V r U(r, t) = uAV(r, t)(43) 
We consider the ansatz 

U(r,t) = W(A(t)-r,t) (44) 
and assume the matrix A(t) to obey the equation 

A T + A T A T = (45) 

This leads to the equation 

^W(e,t)+A-W(^,t)+W(^,t)-A T -V € W(e,t) = uV ( AA T -V^W(U) (46) 
The ansatz 

W(£,t) = $(t)w(£,t) (47) 
where <&(£) has to obey the equation 

+A(t)$(t) = (48) 

leads to 

JUr(£, t) + w(C, t) ■ $ T A T • V 5 w(e, t) = z/V c AA T • V c w(C, t) (49) 
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The matrix $(£), A(t) obeys the equation 



$ + A<l> = 

A + AA = (50) 

We do not pursue the question further, what conditions have to be imposed 
on the matrix G(t) in order to allow for a complete reduction of equation 
(pf9")) to the Burgers equation (possibly in a space with lower dimensions). 
The simplest case is clearly the case of a diagonal matrix 

G(t) = g(t)E (51) 

Then, all matrices $, $ T , A, A T are diagonal matrices and, therefore, coin- 
cide. The matrices 

AA T = $ T A T = \ 2 {t)E (52) 
are proportional to the unit matrix such that a transformation 

dj f^ (<*> 

actually defines w(£, r) as the solution of the multidimensional Burgers equa- 
tion. We mention that a similar treatment has been performed for the case 
of stretched vortices in hydrodynamic flows [Zj 



Summarizing we have considered a forced Burgers equation with a forcing 
linear in the spatial coordinate. We have shown that in the one dimensional 
case this problem can be solved completely. We have discussed the rela- 
tionship to the velocity increment statistics of the Burgers equation with 
translationally invariant forcing. In the multidimensional case we have indi- 
cated the possibility of a similar treatment, which shows an analogy to the 
treatment of stretched vortices in hydrodynamic flows. 
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